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■ In this paper, we consider bosonic reduced Yang-Mills integrals by using some approxima- 

tion schemes, which are a kind of mean field approximation called Gaussian approximation 
and its improved version. We calculate the free energy and the expectation values of various 
operators including Polyakov loop and Wilson loop. Our results nicely match to the exact 
and the numerical results obtained before. Quite good scaling behaviors of the Polyakov 
loop and of the Wilson loop can be seen under the 't Hooft like large N limit for the case 
of the loop length smaller. Then, simple analytic expressions for the loops are obtained. 
Furthermore, we compute the Polyakov loop and the Wilson loop for the case of the loop 
length sufficiently large, where with respect to the Polyakov loop there seems to be no known 
results in appropriate literatures even in numerical calculations. The result of the Wilson 
loop exhibits a strong resemblance to the result simulated for a few smaller values of N in 
the supersymmetric case. 
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1 Introduction 



One of the most exciting topics about nonperturbative aspects of superstring theory or M- 
theory is their various connections to gauge theories. Matrix theory, whose classical action is 
given by that of one- dimensional Yang-Mills theory with maximal supersymmetry, has been 
conjectured as a nonperturbative definition of M-theory by Banks, Fischler, Shenker and 
Susskind (BFSS) [|l|. Also, via its toroidal compactification on a circle and on a two-torus 
in the procedure of Taylor , it leads to a proposal for a nonperturbative definition of type 
IIA superstring theory Q and type IIB superstring theory 0, |], respectively. They are 
given by two- and three-dimensional maximally supersymmetric Yang-Mills (SYM) theories. 
In addition, another definition of type IIB superstring theory, which takes the form of the 
above supersymmetric Yang-Mills theory completely reduced to a point, has been proposed 
by Ishibashi, Kawai, Kitazawa and Tsuchiya (IKKT) 0. 

However, in order to definitely extract nonperturbative aspects of string theory from these 
matrix string models, we need to know also about those of SYM theories. In general, it still 
remains to be an extremely difficult problem, so until now the analysis has been done in the 
limited cases. In energy sufficiently higher than the string scale, there exists the region such 
that physics can be well described by a perturbative expansion around the specific instanton 
configuration of the SYM theory which describes the world sheet geometry in the string 
scattering process p[. In the case of the IIA matrix string theory, some interesting results 
coming from nonperturbative natures of strings have been discovered there — existence of 
the minimal distance, limitation of the number of short strings and pair creation/annihilation 
of D-particles in intermediate states of the scattering process. Also, by mapping the SYM 
theories into cohomological field theories, the partition functions and some special operators, 
which must belong to the BRST-cohomology class in the cohomological theories, can be 
exactly calculated |P, [l^, |Tl|, [l^. There are some approaches to Matrix theory from the 



point of view of a generalized version of the conformal symmetry in four-dimensional SYM 
theory |T^. Besides the analytical computations as above, the numerical analysis has been 



proceeded in IKKT model [0, [T|, [T§, [T^, |g, ^ and in BFSS Matrix theory [^. 

Recently, Matrix theory has been analyzed in the approach based on Gaussian approx- 
imation [^, where it is reported that the Gaussian approximation nicely captures the 
qualitative strong-coupling behavior in various simple models — zero-dimensional 0^-theory, 
M = 2 supersymmetric quantum mechanical systems — and in Matrix theory case it gives 
the results fitting the predictions from the conjectured AdS/CFT duality p^ . 

Here, we concentrate into the bosonic Yang-Mills integrals in zero-dimension as a prepa- 
ration to the analysis of the supersymmetric case (IKKT model), and calculate various 
quantities by using the Gaussian approximation and by its improved version. We calculate 
the free energy, the so-called space-time extent^ and the expectation values of Polyakov loop 
and Wilson loop with a square-shaped contour. Our results match well with the results 
of either exact or numerical calculations reported before. Furthermore, by means of the 



^In this paper, the name "space-time extent" is used for the quantity i / (itr(Xj 



1 



improved mean field approximation, we calculate the expectation values of the Polyakov 
loop and the Wilson loop in the case of the length of the loop sufficiently large, which is 
guessed to reproduce the correct behaviors at least qualitatively. Since in our knowledge the 
expectation values of the loops in this case have not ever been evaluated either analytically 
or numerically in literatures, the approximation scheme presented here might give some new 
insights. 

The paper is organized as follows. In section 2, we explain two approximation schemes 
— Gaussian approximation and improved mean field approximation — by applying them to 
simple (/)^-integral. We calculate the free energy and the correlator {e^^'^). In the Gaussian 
approximation, it can be seen that the first few terms in the expansion around the Gaussian 
classical action yield a good approximation of the exact result. For the correlator (e*^*^), 
this calculation can be trusted in the case of L small, while, as L increases, higher order 
terms becomes more and more relevant. So we have to consider another approximation 
scheme suitable for the case of L large. The Gaussian approximation can be interpreted 
as a kind of mean field approximation. From the point of view, we can improve the mean 
field approximation to be appropriate for treating the large L case. Then we find a precise 
correspondence between the solutions of the mean field equation and of the saddle point 
equation. As the consequence, it turns out that the improved mean field approximation gives 
the result which nicely captures the essential features of the asymptotic behavior evaluated 
by the saddle point method. In section 3, we review briefiy some known results about the 
reduced Yang-Mills integrals for referring them later. Also, an exact result for the space-time 
extent in the SU (2) case is added, which has not ever appeared in the literatures as far as we 
know. Section 4 is devoted to analysis of the reduced Yang-Mills integrals by means of the 
Gaussian approximation. Our results for the free energy and for the space-time extent fit 
well the exact results in the case of the dimensionality of space-time large. This is consistent 
with the fact that the Gaussian approximation is regarded as the mean field approximation. 
Furthermore, with respect to the Polyakov loop and the Wilson loop, our results match 
better with the numerical results given in ref. [jl^ in the region of the length of the loop 
smaller. In the case of general N, the formulas which we have obtained are relatively lengthy. 
However, after taking the 't Hooft like large limit, they become remarkably simplified. 
Looking at the several results obtained for various values of A^ together, we can observe quite 
good scalings in particular in the region A^ > 48. In section 5, we continue the analysis using 
the improved mean field approximation. The result for the space-time extent shows a closer 
approximation to the exact result than that of the Gaussian approximation in the SU{2) 
case. Also, the one-point functions of the Polyakov loop and the Wilson loop are obtained 
in the region of L large. In our knowledge, they have not ever been evaluated either exactly 
or numerically in appropriate literatures. On the other hand, in the supersymmetric case, 
the Wilson loop amplitude for a few smaller values of A^ has been numerically computed up 
to larger L in ref. [1^]. Here we find a strong resemblance between the result there and 



^M. Staudacher informed us that the authors of ref. had simulated the Wilson loop also in the 

bosonic Yang-Mills integrals, and that the similar result as in the supersymmetric case has been obtained 
lESl. It is consistent to our result. We would like to thank M. Staudacher for his kindness. 



2 



ours. We summarize the results obtained here and discuss about possible future directions 
in section 6. In appendix A, in order to confirm the correspondence between the solutions of 
the mean field equation and those of the saddle point equation, which is mentioned in section 
2, we give another evidence which realizes the correspondence by investigating 0^-integral. 
Appendix B is devoted to a detailed derivation of the 't Hooft limit for the Polyakov loop 
and the Wilson loop. 



2 Gaussian and Mean Field Approximations 

In this section, we explain two approximation schemes, which we will apply to Yang-Mills 
integrals later, by using a simple example (0'^-integral in zero-dimension). The first scheme 
is the so-called Gaussian approximation, which is discussed in the case of various supersym- 



metric quantum mechanical systems in ref. |22, 23]. The Gaussian approximation can be 



regarded as a kind of mean field approximation, and then we can consider some improvement 
for the calculation of various correlators. This is the second one. We will call it the improved 
mean field approximation. 



2.1 A Simple Example 

We start with 0^-integral in zero- dimension with the action: S = -^<P'^- The partition 
function Z can be exactly calculated in terms of the Gamma-function: 

-s [9T^f^ 



Now, as a scheme which reproduces this result approximately, we consider the expansion 
around the Gaussian action 5*0 = ^4>'^- 

^^-So^-is-So)^Zo{e-^'-'%, (2.2) 



— oo 



where Zo and (■ • ■)q denote the partition function and the expectation value in the Gaussian 
theory. The width a is to be determined later. The free energy F = — \nZ is given by the 
form of the Cumulant expansion: 

F = Fo-J2 -^{(S - Sor)c,o- (2.3) 

The subscript C means the connected expectation value in the Gaussian theory. The first 
few terms of the expansion are given by 

Fo = -InZo = -^ln(27ra2). 



{S - So)o 



2 

3a^ 1 



4((^-^o)V,o^-3^ + |,-^. (2.4) 
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Of course, the free energy F is independent of a by definition. However, if we truncate 
the series ( |2.3| ) at some order, the value of the truncated series depends on a. Generically 
summing up the whole series (|2.3| ) seems to be a formidable task, thus it will be best to 
choose the value a so that the series exhibits a sufficiently fast convergence into the limit F 
if it's possible and then to evaluate the truncated series at the optimized value a. 
Now let us determine a by the equation: 

^ {Fo + {S-So)o) = 0. (2.5) 



It can be interpreted as the variational method, because the inequality F < Fq + {S — So)o 
holds and we consider the a minimizing the r.h.s. of the inequality^. Eq. ( p.5|) fixes a as 



= then the free energy becomes 

^ = In^? - ^ In ^ - ^ - ^ + 0(((S - So)%,o), (2.6) 

where the third and fourth terms represent the contribution from (5* — 5*0)0 ^-^d — |((S' — 
5'o)^)c,o respectively. Comparing the exact result 

1 rfi) 1 

F = — In ^ - In = — In c/ - 0.94144 ■ ■ ■ , (2.7) 

2 v2 2 



with the result up to the first order in 5 — S'q:"— ^Ingf — 0.89428- ■ ■" and that up to the 



■ " i 

"2 

second order: "— ^ \ng — 0.97761 ■ ■ ■" , it seems that the series ( |2.6| ) tends to converge into the 
exact value. 

Correlator (e*^*^) Next, we consider the expectation value of the operator e*^*^ (which 
corresponds to Polyakov loop or Wilson loop in reduced Yang-Mills integrals we will discuss 
later): 

1 /-oo lp-{S-So)piL(j>\ 



Expanding around the Gaussian theory, we have 



(e^^<^) = (e*^^)o + E ^^((5 - 5o)"e^^^)c7,o. (2.9) 



„=i 



At the value a = y g'^/S satisfying eq. ( p.5|) , the first few terms give 

36 



Eq. (|2.5|) can also be regarded as mean field approximation as we will show in the next subsection. 
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where the second and third terms represent the contributions from the first and second 
order terms in the (5* — 5'o)-expansion. Note that {{S — So)"'e''^'^)co {n > 0) has the form 

Pn{L'^g)e 2^3. p„(a;) is a polynomial of x of the degree 2n with no constant term. We can 
easily see this from a direct calculation except the point that the polynomial has no constant 
term. Including no constant term is understood from the fact that the connected correlator 
{{S — So)"'e'^^'^)c,Q vanishes as L — 0. Thus it can be expected that the expansion ( |2.10| ) 
gives a reasonable result at least qualitatively for the region of L'^g small. 

However, for the case of L^g large, since higher order terms become more and more 
dominant, we can not trust results obtained from the above expansion. In fact, in this case 
we can evaluate the asymptotic behavior of (e*^*^) by using the saddle point method for the 
integral d(j) e~^e'^^'^. The saddle point equation has three solutions: 



5W6(^^^2^|l/3 



(2.11^ 



We deform the integration contour so as to pass the two points (f)± along the steepest descent 
directions and evaluate the Gaussian integrals^ The result after dividing by the partition 
function (E.ll) is 



r(i) 



3 {L^y 



/3 



cos 



(2.12) 



The result ( p. 12 ) exhibits qualitatively distinct behavior from the result obtained from the 
first few terms in the Gaussian approximation (|2.10|) . In the next subsection, we discuss 
some improvement of the approximation which reproduces the behavior ( p.l2| ) at least qual- 
itatively. 



2.2 Improved Mean Field Approximation 

Bearing in mind the improvement, we begin with giving another interpretation of the Gaus- 
sian approximation as mean field approximation. Let us consider the mean field approxima- 
tion to the 0^-integral by replacing 0^ in the 0^-action with the expectation value (0^). The 
mean field action is 

+ + i (2.13) 

where the factor "6" stands for the number of ways of contracting two 0's in the 0^-action, 
and the constant A was introduced for later convenience. The width of the Gaussian action 
a is related to the two-point function as 



-, = -,{<P')- (2.14) 
g'^ 



For discussion in the next subsection, we note that the integrand at the saddle point 
which exhibits the behavior of blowing up as L — )■ oo. 
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The partition function is written as 

Z = 



(2.15) 



The partition function and the expectation values under the mean field theory with the 
classical action Sm are denoted by Zm and (■ ■ ■)m- Since the factor e"*^"^""^*^^ represents the 
difference between the original theory and the mean field theory, we want to take 



-{S-Sm) 



M 



(2.16) 



by choosing the parameters a and A, so that the mean field theory realizes the original 
theory. Also, for the equivalence of both theories[| we require 

(0') = (0')m. (2.17) 

The two conditions ( 2.16| ) and ( 2.17 ) determine the parameters a and A. By combining with 
eq. ( p^.l4D , the condition ( |2.17| ) gives the same value of a as from the variational method 
(^. Also, the condition (|2A6| ) means En=i ^-^{(S - SmT)c,m = 0. It is noted that 1/A 
appears only in the first term of the l.h.s. of this equation because considering the connected 
correlators. From this equation, A is given by 

-T.—^{iS-S',,r)c,M (2.18) 



A 



n=l 



IhZm = — |ln(27ro"^ 



leads 



where we put S*^^ = -^<t>^- The free energy F = —\iaZ = — luzjm 
the identical result with the variational method applied to ( |2.3| ). 

Next we consider some improvement of the above mean field approximation in the case of 
the unnormalized expectation value of the operator O. Let us repeat the same argument by 
regarding the unnormalized expectation value (O)' = J d(j) e~^0 as the partition function of 
a theory with e~^0 being the Boltzmann weight. We take the following mean field action: 

1 ^ , , o 1 , o 1 1 



s 



M 



4g' 



1 

+ - 
A 



2a2 



1 

+ 

A 



(2.19) 

Here, (■ ■ ■)'~' stands for the expectation value under the Boltzmann weight e~^0, i.e. for 
arbitrary operator A, 



(A) 



o 



Jd(f)e-^OA (OA) 



Jd(f)e-30 {O) 
As the result of the same argument as before, we obtain 
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1 
A 



E 

m=l 



m! 



3^2 

■((5-5y™)c,M-ln 



{0)m - 
1 



1 + 



{O) 



E 

M n=l 



-{o{s-s',,r) 



C,M 



[2.20) 

[2.21) 
[2.22) 



where S'j^j = 



^Note that the equivalence stated here is according only to the partition function and the two point 
function, not to all correlators. 
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Correlator (e*^"^) Let us apply this method to the case of O = e^^'^. From the condition 

= -\L^g^ + (2.23) 

where x = L'^a'^. In the case of L^g large, this equation can be solved iteratively and the 
following three solutions are obtained: 

1/3 



xo = -Q^V)'^' + ^ + 0((LV^)-^/^). (2.24) 



Corresponding to each solution of x, A is determined by eq. ( p.22|) . In the form of the 
{S — S'^)-expansion, it is expressed as 



A 12 12 
up to the first order, and 



^ ^ X - 4 + 0((LV^)-^/3) + OiiS - S',,]'^ 



1 OQ 1 Q 



(2.25) 



(2.26) 



up to the second order. We will denote the three A's corresponding to x± and xq by A± and 
Ao, respectively. 

Here, we have three sets of the solutions (5", A). There is a problem — which solution 
should we take? Also, does it make sense to adopt solutions more than one? If it does, how 
should we combine them? In this case, we can propose a prescription based on the following 
observation. The unnormalized correlation function is expressed in terms of x and A as 



— v2vnr exp 
Ij 



1 1 
"2^^" A 



(2.27) 



Examining the exponentiated factor in this equation, we are tempted to claim that these 
solutions reflect the structure of the solutions of the saddle point equation ( |2.11| ). Actually, 
xq and Ao lead to an unphysical amplitude which blows up as exp[c(Ly^)''/^] with c being a 
positive constant when L — oo. This behavior is very close to that of the integrand at the 
saddle point 0o (See footnote ^ in the previous subsection.). Furthermore, we can see that 
both of (x+,A+) and 0_ give the similar behavior as exp[— e'^*'^^c(Ly^)^/^] and that both 
of (a:_,A_) and 0+ lead exp[— e~'^*/^c(Ly^)''/^]. In particular, the solutions (x. A) exactly 
reproduce the phase factors and the power of L^^g in the saddle point values of —S + iL(j). 
From this fact, it will be plausible to assume that each solution (x. A) certainly corresponds 
to each saddle point solution, although at present we have not found out the definite reason 
for the correspondence^. If accepting this assumption, we can proceed further. As mentioned 

^In order to confirm this assumption further, we give another evidence for the case of 0^-integral in 
appendix A. 
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above, because the solution (xq, Aq) lead the unphysical solution, we discard it. Also, since 
the unnormalized expectation value (e*-^*^)' = (cos(L0))' is real, we combine the two solutions 
[x±, X±) with an same weight. Here, we determine the weight from the assumption. Let us 
take the combination same as what appears in the saddle point method. That is, we simply 
sum up the contributions from the two solutions with the weight one: 



{e'^^y = -^2t:x+ exp 



1 

X J 

2 ^ 



A, 



+ —d2TiX- exp 



1 

IT 



(2.28) 



Finally, dividing by the partition function, we obtain the expression of the correlator. 
Then, for the precise cancellation of the vacuum graphs between the unnormalized correlator 
and the partition function, we need to use the result of the partition function up to the order 
same as that of the unnormalized correlator^. Thus, for the result up to the first (second) 
order we use the result of the free energy up to the first (second) order. The final expression 
is 



x[l + 0((LV^)-^/3)+0((5 



COS 



4/3 



J A 



\n+l 



with up to the first order (ra = 1): Ui 

29 „, _ 41 



36' 



M) 

n 

24' 



and up to the second order 



n 



;2.29) 
= 2): 



U2 



27' 



V2 



72- 



Now let us compare these with the result of the saddle point method (|2.12| ). First, the 
power behavior l/(Ly^)^/^ and the power in the exponential and the cosine (Ly^)^/^ just 
coincide. The coefficient of the exponential decay is 

il . 3-1/3 = 0.31779 ■ ■ ■ up to the 1st order 

41 0-1/3 

F 



0.39483 



(2.30) 



0.375 



up to the 2nd order 
saddle point method, 

where the second order result approaches closer to the saddle point result than the first 
order result. Also, the coefficient appearing in the argument of the cosine exhibits the 
similar behavior. On the other hand, the constant factor in front of the whole expression 
does not show a good result as long as looking at the first two orders: 

2 ■ 3i/i2e-i 



2 ■ 3i/i2e-i 



1.3668 ■• 
0.74873 ■ 
1.1289- •• 



up to the 1st order 
up to the 2nd order 
saddle point method. 



:2.3r 



We need further examination of higher orders for convergence of the constant factor. The 
decay coefficient is determined by the first term in (|2.25| ) (or (p.26|) ) which is the leading in 
the case large. On the other hand, the constant factor is by the second term which 
is the sub leading. From this point, we can understand that the decay coefficient converges 
faster than the constant factor. We can conclude that our scheme quite nicely reproduces 
the qualitative behavior of the saddle point result. 

^Namely, it guarantees (1) = 1 in the result up to every order. This consideration becomes more relevant 
when treating the system with more degrees of freedom. 
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3 Reduced Yang-Mills Integrals — Exact Results 



In this section, we give some explanations about reduced Yang-Mills integrals, before ap- 
plying the method discussed in the previous section. We review some known results as well 
as add an exact result for the space-time extent which have not been seen in appropriate 
references. 

We consider the following bosonic Yang-Mills integral: 



D 

-S{X) 



nMX,]e-^(^), (3.1) 



with the Euclidean classical action: S{X) = — ^tr[X^, Xj,] [X^, X^]. The variables X^'s 
are traceless hermitian matrices of the size N x N. The indices yU and u run from 1 to D. 
The normalization of the measure is determined by|^ /[c^-^] e"*''"^ = 1. At first sight, this 
integral (|3.1|) seems to lead an ill-defined result due to the integration over the infinite range 



along the flat directions. However, as shown in ref. u3, 16], in the N = 2 case the integral 



can be performed exactly and it turns out to give the finite result when D > A: 

[ oo {D <A) 

Zd,n=2=\ i.,fD r(f)r(f-i)r(f-i) , . (3.2) 

As is understood from intermediate steps in the integral, the finiteness is thanks to the 
suppression factor generated after integrating out the other variables than the variables 
along the fiat directions. Furthermore, along the similar line we find the exact expression 
with respect to the square of the so-called space-time extent (|tr(X^)^): 



;^tr(X,)^) = |(D-2)^A|_^ (3.3) 



for D > 6. It is not well-defined for the case oi D < 6. This result had not been obtained in 
literatures as far as we know. 

In the case of X > 2, there are some insights from perturbative analysis. Let us start 
with considering the polar decomposition X^ = V^A^V^, where A^'s are matrices of the 
eigenvalues of X^'s: A^ = diag(A^, ■ ■ ■ , A^^), and V^'s are unitary matrices. The integrals 
of the unitary matrices are performed perturbatively by expanding around the unit matrix. 



As shown in ref. p6[, from the formula of Zjj ^ after integrating V^'s at the one-loop level, 
powercounting with respect to the A^-integrals leads to the condition of the convergence of 
the integral for the large separation among the eigenvalues: 

X > — (3.4) 
D-2 ^ ' 



^"^This is nothing but the same normahzation of the measure in ref. p5[ . It can be easily seen by expanding 
X by the basis T-^'s normalized as tiT^T^ = ^5ab- X = 'ZaI'i^ X^T^. 
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For the eigenvalue density of one of the D-matrices (say, X/)): p(A) = {jj J2i=i ^{^ ~ ^d))^ ^ 
similar but more careful consideration about the A^-integrals leads to the following asymp- 
totic behavior: 

p(A) ~ A-2^(^-2)+3^-^ (3.5) 



51) are consistent with the 



which has been derived in ref. [T6|. The formulas (|3.4| ) and 
results (13) and in the = 2 case. 

In the bosonic Yang- Mills integrals, we can also analyze by the 1/D-expansion method 
as in ref. l2^^^. For a later discussion, we show the result for the square of the space-time 
extent: 



(^tr(X,)^) 



iDNg^ - 1 



1 + 



1 

D 



1 



0{D- 



(3.6) 



4 Reduced Yang-Mills Integrals — Gaussian Approx- 
imation 

In this section, we apply the Gaussian approximation explained in section 2.1 to the reduced 
Yang- Mills integrals. With respect to the partition function and to the space-time extent, 
our results nicely fit the known exact results of the N = 2 case when D is large. Also, 
we compute the one-point functions for the Polyakov loop and for the Wilson loop with a 
square-shaped contour. Our results turn out to reproduce well the numerical results in ref. 



17] when the length of the loops is smaller. 



4.1 Partition Function and Space-Time Extent 

First, we start with the partition function. Expanding around the Gaussian action: 5*0 = 
^tr(X^)^, we get the same expression as ( |2.3D with 

Fo = -^{N^-l)D\na\ 

{S - So)o = f-,NiN' - l)D{D -1)-\{N'- 1)D, 

-\{{S - S,f)c,, = -^,N\N' - 1)D{D - 1){AD - 1) 

^^The 1/D-expansion method enables to evaluate various correlators and their large- TV scaling properties 
systematically as presented in ref. At present, unfortunately it seems hard to systematically calculate 

highly complicated composite operators such as Polyakov loop and Wilson loop, which we discuss in this 
paper, by using the l/D-expansion formalism. The Gaussian and improved mean field approximations 
are not systematic from the point of view of the l/I?-expansion, which is a weak point of these methods. 
However, they enable to calculate those loop amplitudes and give a closed form at each order with respect to 
{S — S'o)-expansion and (5 — S'^^)-expansion. As we will see, the results show good agreement with numerical 
results. 
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4 1 

^ N{N^-1)D{D-1)--{N^-1)D. [A.l] 



4^2- V- -/V 4 
As the result of the variational method (|2.5|) , a is determined as 



2^^ 



As we saw in section 2, this approximation can be regarded as a kind of mean field approx- 
imation. So it is interesting to compare this with the exact result in the N = 2 case ( |3.2| ) 
when D is large. The exact result behaves as 

- In Zd,n=2 = ^D\nD - ^D{\ng^ + 1) - 2\n2 + O (^^y (4.3) 

while the result by the Gaussian approximation becomes correspondingly 

3 3 91 / 1 \ 

F=^DlnD- ^D{\ng' + 1) - _ + O (- j + 0{{{S - So)')c,o). (4.4) 

The first two terms in both formulas completely coincide. Also the 0(D°)-terms are —2 In 2 = 
—1.3862 ■ ■ ■ and — = —1.3125 ■ ■ ■. It seems that the Gaussian approximation quite nicely 
reproduces the exact result in the large D case. Furthermore, we can compare with the 
numerical result for smaller D and A^'s (Table 1 in the second paper of ref. ||16|), from which 
the values of -F(D,Af) = — InZ + ^D{N'^ — 1) ln(2(yf^) are read off as 

F(3,4)=4.98---, F(3,5) = 13.8---, F(3,6) = 26.1 ■ ■ ■ , 

F(4,3) = 6.27 ■ ■ ■ , F(4,4) = 17.4 ■ ■ ■ , F(4,5) = 33.8 ■ ■ ■ , F(4,6) = 56.2 ■ ■ • . ^ ' ' 

Let us look the quantity '^{d,n) = (-^(D,Af) ~ ^(d,n))/ F'(d,n) where FI^d,n) the quantity in 
the Gaussian approximation corresponding to F(^£),n)- Then, 

A(3,4) = 0.59 ■ ■ ■ , A(3,5) =0.21---, A(3,6) =0.11---, 

A(4,3) =0.21---, A(4,4) = 0.063 ■■• , A(4,5) = 0.036 A(4,6) = 0.021 ••■ . 

This shows a tendency of better agreement not only for larger D but also for larger A^. 
Next, let us examine the square of the space-time extent (-^tr(X^)^). We have 

T-Mx.f) = {-Mx,r)o + E - sor-Mx,r)c,o, (4.7) 



with 



-{{S - 5o)^tr(X,)2)c,o = -^(A^^ - 1)D{D - 1) + y 

1 1 fT^O 

-{{S - Sor-tT{X,r)c,o = Y^,N{N' - l)D{D - 1)(4D - 1) 

'^\N^-l)DiD-l) + ^-^D. (A, 



4^2^ ' ' ' 2 N 



11 



At the value ([4 .21), the first order term vanishes, and we obtain the expression 



(]^tr(X,)^) 



Ng^ - 1 

^ 2{D-i) m 



D 



2{D 



+ OiiS~Sof) 



(4.9) 



When D is large, the second order term is suppressed by the factor comparing to the 
zeroth order term. It is consistent to the picture of the Gaussian approximation as the mean 
field approximation. In the N = 2 case, this has the following large D expansion: 



:^tr(X,)^) = '-IVD 



2 
D 



l + ^ + OiD~^) + Oi{S-SoY) 



On the other hand, the exact result (|3.3|) does 

(^tr(X,)^) = ^v^ 



1 + ^ + 0(1.-) 



(4.10) 



(4.11) 



We find some difference between the 0{D ^/^)-terms. As will be seen in the next section, 
the situation becomes better when applying the improved mean field approximation. 



4.2 Polyakov Loop 

Here we consider the expectation value of the operator of a loop of the length L winding in 
one direction (say, the first direction) P{L) = -^tre*^^^, which we will call Polyakov loop. 
The expectation value is expanded as 

(P(L)) = (P(L))o + E ^^((^ - SorP{L))c,o. (4.12) 

n=l ^• 

For calculating each term in the r.h.s. it convenient to consider the Gaussian integral over the 
hermitian matrices (including the trace part) and to use the orthogonal polynomial method 
for the integrals with respect to the eigenvalues. Let us introduce the hermitian matrix Y 
by adding the trace part to the traceless hermitian matrix X: Y = X + V'^^'^n, where 

y eH. Also, the measure [dY] is normalized by /[ciF] e"*''^^ = 1. We consider the following 
expectation value in the hermitian Gaussian integral: 

The Gaussian weight can be factorized into the product of the trace part and the traceless 
part. After integrating out the trace part, eq. ( [4.13| ) turns out to be related to the Gaussian 
expectation value of the Polyakov loop as 

(ltre^^^)o = e-'i^(P(L))o. (4.14) 
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Also, the l.h.s. of this equation reduces to the integrals with respect to the eigenvalues 
{Xi} of Y, and it can be easily evaluated by using the orthogonal polynomial method [P^ p|. 
It is translated into the language in the quantum mechanical system of a harmonic oscilla- 
tor. Connected correlators among the [/(A^) -invariant operators tTFk{Y) (fc = 1,2, ■ ■ ■) are 
expressed as 



N-l 



{tTF,{Y)),= J2{n\F,iX)\n), 



n=0 



(trFi(F) trF2(y))c,o = Tr[njvFi(A)(l - Un)F2{\)], 

(trFi(r) trF2(r) trF3(r))c,o = Tt[UnFi{X)F2{X)Fs{X) -U^FiCmNF2{X)Fs{X) 
-Fi(A)n^F2(A)n^F3(A) - n^Fi(A)F2(A)nA,F3(A) 

+n^Fi(A)n^F2(A)n^F3(A) + n^Fi(A)n^F3(A)n^F2(A)], 

(4.15) 

where the creation and annihilation operators a and appearing in A as A = -^(o + S^) 
satisfy [a, a"'"] = 1. The states {|n)}„=o,i,--- form an orthonormal basis in the Fock space of 
the system of the quantum harmonic oscillator: d\n) = \n — l)^/n, d^\n) = \n + l)^/n + 1. 
"Tr" means the trace operation over the infinite dimensional Fock space. Also, "11 at" stands 
for the projection operator into the A^-dimensional space: IItv = J2n=o -^y using the 



first formula in eqs. ( [4.15|) we find 

, 1 



tre^^O 



-F 1~N,2; 



(4.16) 



where F is the confluent hypergeometric function: 

a{a + 1) ■ 



F{a,^-z) = Y: 



ia + n- l) z''' 



;^o/?(/3 + l)---(/3 + n-l) n! 
Thus from eqs. ( |4.14| ) and ( ^.16| ), we obtain 



(^(^))c 



AT, 2; 



(4.17) 



First Order Term Next, we compute the first order term, namely — ((S* — So)P{L))c,o- 



Let us consider {SoP{L))o: (S'oP(L))c 



D-l 
D 



{So)o{P{L))o + (^tr(Xi)2p(L))o. We pass to 



the F-integral in order to evaluate (tr(Xi)^P(L))o. Then, 



r2 2 1 1 1 



{SoPiL))o = e—{-tTY'-tTe'^')c,o + {^tTY%{PiL)) 



+ 



D - 1 



D 



N 

{So)o{P{L))o 



(7^ 



2 



(4.18) 



^^Unfortunately, in the standard framework of the l/D-expansion m ref. 
orthogonal polynomial method, which is the difficulty pointed out before. 



I, it seems hard to apply the 
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By making use of the second formula in eqs. (|4.15| ), the connected correlator in the first 
term can be calculated. Evaluating {SP{L))q similarly, we eventually obtain 

^4 \ 



-{(S-So)PiL))c,o 



N{D - 1) - 1 



{N - 1)F h - N, 3; ^ j +{N + 1)F il - N, 3; -^j 



AN 



(4.19) 



Second Order Term For the second order term ^{{S — Sq)^P{L))c\o, we can also compute 
in the similar manner as in the first order case. After a relatively long but straightforward 
calculation, we arrive at the following result: 

l{{S-SorP{L))c,o = 



g 4 N ' 



1 - -N(D - 1) - _(JV - 3)(B - 1) + —^N'(D - If 



4 



1 - ^-N{D - 1) + ^(iV + 3)(D - 1) + ^^N\D - If 



32N 



32N 



4^4 



r 4 

+ e ^(i-iv) — -(L)-l) 



N + l\ ( L^a^' 
3 P 2-iV,4;^ 



— e 



2 2 1 

'-^il-—') 

4 JV'' 



3a^ 



F 3 - iV, 4: 



3 P^-^'^^^ 



0" 



iV(D - 1) + —N{2N - 1){D - 1) + —N\D - 1)(4D - 3) 

^ 2 



l--iV(D-l) 



— e 



4Ar 

'2^2 r 



16iV2 



2 P 2-^'^^^ 



3a' 



a 



a 



N{D - 1) + — A^(2A^ + 1){D - 1) + ^A^l/^ - 1)(4D - 3) 



x_.+ l-^iV(D-l)' 



AN 



16iV2 



iV + 1 
2 



+ 



3g' 



F 1~N,3; 



1 - ^NiD - 1) + -N^DiD - 1)| ^ + |1 - ^iVP - 1) 



32iV2 
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x(P(L))o. 



(4.20) 



Gaussian Approximation Let us evaluate the first three terms in the expansion (^4.12|) , 
i.e. eqs. { ^1% , i ^J% and (|]20|), at the value (|]2]). The first order term ( ^19[) vanishes, 
and the second order term ([4.20|) takes the form 



X 



+ 





CO = e" 


j N -3 


LV^ " 


\N^{D - 1) 


32A^ 




LV^ ' 


' N^D - 1] 


32N 


1 





F ( 3 - A^, 5; 



F I 




2N{D - 1) 4A^ V 2 

3N + 1 L / N + 1 
'2N{D~1) AN \ 2 
2 L^a^ 




N,5; 



F ( 2 - iV, 4; 



1 4iV 



(4.21) 



where 



29^ 



We can observe that the factor exists in every term in the second 



order result. Again, it is consistent with our picture as the mean field approximation. 

In the case of D = 4, there are some numerical results about the Polyakov loop and the 
Wilson loop reported in ref. 
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So we can compare our result with the numerical one. Let 
us take Ng'^ = 48 and consider the quantity for various values of N keeping Ng'^ fixed, which 
gives the same setting as in ref. [0. Figure 13 in ref. JT^ is the result to be compared 
with ours. The variable k/y/g appearing there corresponds to L in our setting. In Fig. 1 
we show the result of the zeroth order alone and that summed up to the second order for 
the Polyakov loop when N = 16. It can be seen that our result up to the second order 
nicely reproduces the numerical result for the region of smaller L (up to about 1.0). There 
are some differences between them in larger L, where in our analysis higher order terms are 
considered to become more and more important. (See also Fig. 2.) 

T7] exhibits a really good scaling behavior against various 
First let us consider the limit A — oo with Ng"^ 
fixed. In this limit, it turns out that the confiuent geometric functions reduce to the Bessel 



Furthermore, the result in ref. 
A's. It can be seen also in our result. 
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Figure 1: Polyakov loop in the case D = 4 and = 16. The dotted hne shows the 
contribution from the zeroth order alone, and the solid line shows that up to the second 
order. 



functions and that the formula becomes considerably simplified as0. 



{PiL)) 



J4(2v^) - 3 J2(2v^)l + 0{{S - So) 



(4.22) 



where a 



and G stands for the fixed value \Ng. The first and the second 



terms in this formula come from the contributions of the zeroth and the second order terms, 
respectively. In D = 4, we plot this result with = 48 in Fig. 2. Comparing it with the 
curve in Fig. 1, it seems to suggest a good scaling behavior. In fact, we plot our results for 
various values of N in Fig. 3 (A^ = 16, 48, 100, 400 and oo). In Fig. 3 we can see some 
slight differences between the curve of = 16 and the others for L > 3. Among the curves 
for A^ > 48, however the scaling is too good to observe any difference for various A^ in the 
region L < Q. 

Since we have obtained the concrete formula for general D and A^, investigating the 
behavior for other D's can be immediately done. For example, we plot the D = 10 case 
and the D = 26 case for various A^'s in Fig. 4. Both casess still exhibit quite good scaling 
behaviors. As D increases, the curve tends to broaden and to make the amplitude of the 
oscillation decreased. We can understand the tendency by seeing the large A^ result ( ^4.22|) . 
The former can be understood from the fact that L is contained only through the combination 
of a, and the latter from the suppression of the second term when D being large. 

Thus we have seen that the Gaussian approximation reproduces well the numerical result 
of the Polyakov loop in smaller L region. How about the Polyakov loop in larger L region? 
From the same reason as the case of the simple example in section 2, the Gaussian approxi- 
mation can not be trusted. In section 5, we will consider it by using the improved mean field 
approximation, and derive the asymptotic behavior of the Polyakov loop approximately. 
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We give some detailed explanation with respect to the derivation of this limit in appendix B. 
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Figure 2: Polyakov loop in the case D = A and = oo. The dotted hne shows the 
contribution from the zeroth order alone, and the solid line shows that up to the second 
order. The gray boxes indicate some points of the numerical data for the = 48 case in 
Figure 13 in ref. |1^ . 




Figure 3: Polyakov loops up to the second order for {D,N) = (4,16), (4,48), (4,100), 
(4,400) and (4, oo). The gray line shows the case (4, 16), the dotted line (4,48), the dashed 
line (4, 100), the dot-dashed line (4,400), and the solid line (4, oo). 
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Figure 4: Polyakov loops up to the second order for {D,N) = (10,48), (10, 100), (10,400), 
(10, oo), (26,48), (26,100), (26,400) and (26, oo). The dotted hne shows the case (10,48), 
the dashed hne (10, 100), the dot-dashed hne (10, 400), the sohd hne (10, oo), the gray dotted 
hne (26,48), the gray dashed hne (26, 100), the gray dot-dashed hne (26,400), and the gray 
sohd hne (26, oo). 

4.3 Wilson Loop 

We consider the expectation value of the rectangular Wilson loop of the size Lx L: W{L) = 
1.1y(^q^lxi ^iLX2 ^-iLXi ^-iLX2y Here we evaluate the first two terms in the expansion: 

{W{L)) = {W{L))o + E - SorW{L))c,o. (4.23) 

Leading Order Term First let us evaluate the leading (zeroth) order term {W{L))q, 
which is written as 

(#(L))o = j^IiL),,kiIiL)jku, (4.24) 

with I{L)ijki = {{e'^^^^)ij{e~'^^^^)ki)o- From the invariance of the measure [dXi] under Xi — > 
UXiW with U being an arbitrary SU{N) matrix, the quantity I[L)ijki must satisfy 

I{L),^ki = Uu'{U^)y,Ukk'{U%iI{L),^,k'v. (4.25) 
It determines the structure of the indices of I{L)ijki as[^ 

I{L\,u = A{L)6,,6u + B{L)6u6ik. (4.26) 

where A{L) and B{L) are S'?7(A^)-invariant. Noting I{L)ijji = and I{L)ijki = 6jk, A{L) 
and B{L) are determined. Thus we can rewrite {W{L))q in terms of correlators in the 

^''in the N — 2 case, the tensor seems to need to be taken into account. However, we do not have to 
worry about it because of the identity: eik^ji — SijSki — SuSjk- 
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Gaussian one-matrix model: 



{W{L))o = 1 - 



iV2 



tre 



iLXi 



N 



tre 



-iLXi \ 



- 1 



(4.27) 



Considering the integral over the hermitian matrix Y as before, the above correlator can be 
evaluated, and we obtain the expression 



{W{L)), 



1 - 



iV2 



m-1 

where f{x) is a polynomial of x given by 

N-l n-1 



1 9 r2„2 



(4.28) 



1 



n=l m=0 

n—m ( ™ 



E 

^=0 



,5(fc!)2(n-fc)! 
1 



E 



-(/!)2(m-/)! 



A;!(n — m + /c)!(m — fc)! \ 2 



(4.29) 



First Order Term Next, we compute the first order term —{{S — Sq)W{L))c,q- After 
integrating out the variables other than Xi and X2, we have 



-{{S-So)W{L))c,o^{N''-l) 

'2 1 



Ag- 



;N{2D - 3) - 1 



{W{L))o 



a 



{{triX,yWiL))o + {tr{X2rWiL))o} 



-^{{XiX,)ij{X2X2)jiW{L))o + ^{{X,)ij{Xi)ki{X2)jk{X2)iiW{L))o. (4.30) 

Noting the S'?7(A^)-transformation property as in the case of I{L)i.jki. the calculation of the 
correlators in the r.h.s. is reduced to the following more fundamental quantities: 

m ^ (^tr((XO^e^^^^)ltre-^^^Oo, 
h{L) ^ (Itr(XO^^tre^^^^ltre-^^^)o, 



h{L) ^ (ltr(Xie'^^^)ltr(Xie-^^^^))o. 



(4.31) 



The result is 



(tr(A-i)'»'(L))„ = (tr(Xj)'»'(L))„ = -a'(N' - 1 



1 9 r2„2 



(4.32) 



19 



1 /V^ — 1 

-(/l(L)-^(L))^ (4.33) 



(A^2 _ i)(Ar2 _4) 



Q /\r2 _ 1 
((Xi),,(Xi)«(X2),-fe(X2)HW^(L))o = -a^ 



4 Ar(Ar2 _ 9) 

"^'^ (7V2_4)(Ar2_9)^^(^) - (Ar2-4)(7V2_9)^^(^) " ^'^ A^^^^^^ 

-(Ar2_i)(7V2_4)(Ar2_9)^^^^)^^^^) 
2 16A^(2A^2 _ 3) 





- l)(A^2 _4)(7V2 


-9) 




Ar3(Ar2 + 6) 






-1)(A^2_4)(^2 


-9) 




20N^ 






-1){N^-4){N^ 


-9) 



(/^(i)) +(jv2_i)(jv2_4)(JV2_9)-''iWA(i) 

+ (J-^(i^2'!9) (/»(^))^- (4-34) 

Now, what is needed for giving the expression of the first order term is to know the three 
quantities fi{L), f2{L) and fz{L)- We pass to the Gaussian integral of the hermitian matrix 
Y to evaluate them. Going along the same line as before, we obtain the following formulas: 



2A^2 . ^ n V 2 iV / V ' ' 2 / 

(n^m) n,m=0 ^ \ / \ / 

-4nF f -n + 1, 1; -^j - '^n^^F i-n + 2, 3; + AnF i-n + 2, 2; ^ 

-2..(-..l,2;^)).(-.,l;^) 

^i^o^H{n-m)\f\[ 2 ) ^ ' 2 j 

2n + 2m + 1 - — F -m, n - m + 1; 



V2;v NJ \ ' 2y 

+ F -m + l,n-m + 2; 

\ 2 J n-m+1 \ ' 2y 



+ ^ — - — J (n — m)(n — m — 1)F j — m, n — m + 1; 

xF -m,n-m + l;— — , (4.35) 
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-n, 1; 
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r2 2 

+ 2(n- 1)— -F I -n + 2,2; 



n^—F -n + 2,3;^— 



1 j F ( -n + 1, 1; ^ j - 4n^F ( -n + 1, 2; ^ 



F -n, 1; - 2nF l-n+l, 2; + 2nF 



N-l 

E 



2^2' 



-n, 1; 



Af~l n-l I 

2EE 



' r 9 9 \ n— m+1 



+ 



+ 



(n + m)F — m, n — m + 1; 



F I — m, n — m + 1] 



2^2' 



2m 



n — 171+ 1 



F —m + 1, n — m + 2; 



2^2' 



{n — m){n — m — 1)F — m, n — m + 1] 



xF —m, n — m + 1; 



+ 



a 



X 



N-l n-l 



e 2 



-2EE 



n=2 m= 



\m\ {n — m + l)!(n — m — 1)\ \ 2 



n — m F —m, n — m; + (n — m)F —m — l,n — m; 



xF — m + 1, n — m + 2; 



LV2 



^2 2 2 -^^-3 



+ ^e-^2E(n + 2)(n+l)^^F(-n-l,2;^^)F 



-n,2;^ , 
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^ ^ n=l m=0 



EE 



■A^ <! F f -n, 1; ) + 2nF ( -n + 1, 2; 
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-n, 1; 



2^2' 



F -mA;—- 



X 




a 



-m, n — m + 1: 



(4.37) 



We used some recursion relations among confluent hypergeometric functions for later conve- 
nience in considering the large N limit. 



Gaussian Approximation We consider the above result at the value ( |4.2| ) to obtain 
the Wilson loop amplitude up to the first order in the Gaussian approximation. Since the 
formula is quite lengthy, we do not write down here directly. Instead, we plot the result in 
the case of various values of D and N below. As in the case of the Polyakov loop, the large 
N limit also yields the result remarkably simplified as(3 



(WiL)) 



1 



a 



[9{U2^))\U2^)f + ^{U2^)Y) + 0{{S - So)'). (4.38) 



The first line and the first term in the second line in the r.h.s. represent the contributions 
from the leading order term and from the first order term, respectively. The factor 



D-l 



m 



the second line means the validity of the interpretation as the mean field approximation. 

In the D = 4 case, it is interesting to compare our result with the numerical calculation 
reported in ref. |T^. Figure 11 in ref. is the result to be compared with ours. As in the 
case of the Polyakov loop, we take Ng"^ = 48 for the common setting. In Fig. 5 we show the 
result of the leading order alone and that summed up to the first order when N = 16. It can 
be seen that the result up to the first order match well with the numerical result. (See also 
Fig. 6.) We similarly plot the result for the case = oo in Fig. 6. Comparing this with 
Fig. 5, we expect a good scaling behavior in the region of smaller L. In fact. Fig. 7 shows 
the behaviors for = 16, 48 and oo together, from which the scaling for smaller L can be 
observed. In particular, the scaling between the cases of = 48 and A = cx) is too good to 
distinguish each curve in the region L < 4. These scaling behaviors are in conformity with 
the numerical result in ref. |l^ . 

As in the case of the Polyakov loop, we depict the behaviors for D = 10 and 26 in Fig. 
8. We can still see quite good scaling behaviors. Also, from Fig. 7 and Fig. 8, it can be seen 
that the curve tends to become broader as D increases. The tendency can be understood 
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See appendix B for detailed explanation. 
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Figure 5: Wilson loop in the case D = 4 and N = 16. The dotted line shows the contribution 
from the zeroth order alone, and the solid line shows that up to the first order. 




Figure 6: Wilson loop for D = 4, iV = oo case. The dotted line shows the contribution from 
the zeroth order alone, and the solid line shows that up to the first order. The gray boxes 
indicate some points of the numerical data for the N = 48 case in Figure 11 in ref. [17 
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Figure 7: Wilson loops up to the first order for {D,N) = (4,16), (4,48) and (4, oo). The 
dotted line shows the case (4, 16), the dashed line (4,48) and the solid line (4, oo). 



5 Reduced Yang-Mills Integrals — Improved Mean Field 
Approximation 

In this section we evaluate the expectation values for the space-time extent, the Polyakov 
loop and the Wilson loop by the method of the improved mean field approximation. For the 
space-time extent, it turns out to give a result closer to the exact value than the Gaussian 
approximation in the case of = 2. Furthermore, for the Polyakov loop and the Wilson 
loop, it yields the asymptotic behaviors for the large loop. Although we consider only the 
first few terms in the expansion around the mean field configuration, as we have done in 
section 2.2, it can be guessed to give the results which reproduce the correct behaviors at 
least qualitatively. 

Let us consider the improved mean field approximation in the reduced Yang-Mills inte- 
grals for the operator O which is U (A^)-invariant and isotropic in D-dimensional space. We 
start with the mean field action 



from the formula in the large N limit ( [4. 381 ), where L appears only through the combination 



of a. 




XI 2 [-{{X^l)^j{X^l)kl)'^{Xu)jk{Xl,)li + {{Xf,)ij{X^)jk)^ {X^)ki{Xy) 



H 



+ 



A 



1 



(5.1) 



Since O is [/(A^)-invariant 



{iXf,)ij{Xf,)M) 



O — 



JjdX) e-^{X^),,{X,)kiO 
JidX) e-^O 



oc 6ii6jk - —6ij6ki. 



(5.2) 
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Figure 8: Wilson loops up to the first order for (D, N) = (10, 16), (10, 48), (10, oo), (26, 16), 
(26,48) and (26, cxd). The dotted line shows the case (10,16), the dashed line (10,48), the 
solid line (10, oo), the gray dotted line (26, 16), the gray dashed line (26,48), and the gray 
solid line (26, oo). 

Also, due to the isotropic nature, the proportional constant in the above is independent of 
IJ,. So the action Sm can be written as the standard form Sm = ^tr(X^)2 + i with 



1 D-1 



^2 g2 m-l'N 

The self-consistency condition (j^ti^X^Y)'^ = {-^ti{X^Y)'^j reads 



1 2 ^ iV^-l {j^tT{X,fO) 



a^^ D-1 m 



M 



{O) 



M 



(5.3) 



(5.4) 



which determines a. As before, A is given by the formula of the same form as ( |2.22| ) with 
S'm = ^tr(X^)2. 



5.1 Space-Time Extent 

Now we consider the case O = -^tr(X^)2. From (|5.4|) , we have 

2 



\ N{D - 1) 



1 + 



1/2 



Up to the second order with respect to S — S'j^.j, X becomes 
^ = N( - DDi I) - 1) - ^ 

3a^ 



A 



^^NiN^ - l)DiD - 1) - liN^ - l)D + |JiV(D - 1) - ^ 



N\N^ - 1)D{D - 1){AD - 1) 



(5.5) 



N'D{D-l) + 0{{S-S'My).{5.Q) 
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Following the steps in section 2.2, we obtain the expectation value 



X exp 



iV^2 ^2 _ ^ 



\2{D-l] 

1 3 

2 ^ A{D - 1; 

1 + 



Ar2 



+ 
2 



3 



(iV2 - 1)D 
n -2" 



+ 



6 



4(D - 1) (A^2 _ 1)/) (iV2_ 1)2/^2 



{N^ - l)D_ 

+ 0{{s - s'm?)]. 

Then, in the case of = 2 and D large, it has the following expansion: 

13 



1 



12D 



+ o{D-') + o{{s - s'^,y 



(5.7) 



(5.8) 



Here, it is noted that this result fits the exact result ( [4.11|) better than the case of the 
Gaussian approximation ( |4.10|) . So, this scheme can be expected to give an improved result 
also for other N. Moreover, there is another comment. If we stopped the computation up 
to the first order, we would end up with the expansion 



1 



tr(X,)2) 



3^? 



I + ^ + 0{D-') + 0{{S - S',,) 



(5.9) 



Thus, we can see that in the SU{2) case the result approaches closer to the exact result as 
increasing the order of the approximation. For general A^, we can compare with the 1/D- 
expansion result ( p.6|) . Normalizing the leading term to the unity, we have the following 
coefficients in the 0(D~^/^)-terms: 



1 



2{Ar2-l) 

2 + 2(Af2-l) 

I J ^ 

6 ^ 



Gaussian approximation up to the 2nd order 
Improved mean field approximation up to the 2nd order 
Improved mean field approximation up to the 1st order 
1/D-expansion method. 



(5.10) 



which exhibits a similar tendency as seen in the N = 2 case. Thus it can be expected that 
the expansion with respect to — S'j^ based on the improved mean field approximation does 
work for general A^. 



5.2 Polyakov Loop 

Here we consider the case O = -^^2^=1 ;^tr(e*^"^^ + e^*^"^''), whose expectation value is 
equal to {P{L)) if 0(D)-symmetry is not broken spontaneously. In the case that N and 
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D are finite, it always holds. Let us consider this case0. For this O, the method in the 
beginning of this section can be used. Now, eq. (|5.4|) leads 



5^^ D-1 ^ 4 i? _ AT, 2; ^ 

We are interested in the case of L^g and x = L'^a'^ large. In this case, the ratio between the 
confluent hypergeometric functions is expanded as 



F(l-iV,3; 



F(l-Ar,2;f) N+1 
so the above equation has the form 



l-(iV-l)^-iV(iV-l)-i + 0(x-^) 



g'^L^ D 1 _ 1 1 1 2N 1 2 

x2 D-1N~ 4N{N + lf~2 WTl'JiTTx^ ^ ^ 

We iteratively solve this equation and get the three solutions: 

x± = e±-^/3(^-^(iV + l)LV) +-NiiN+l)D + 2) + 0iiL^)-'/'), 

^—^{N + 1)LV) + ^N{{N + 1)D + 2) + 0{{L^r'/'). (5.13) 



Xo 



Corresponding to each x, A is determined by 

i = (iV - 1) {^a; - 1 + 0{x'') + 0{{S - S',,)')] (5.14) 
up to the first order with respect to 5* — S'^j, and 

^ ^N-l)\^x~-{N+l){N-9)D--{N + lf ^ ^ 



A 12N 8' ' 8' ' D-\ 2 

+0(x-^)+0((5-5;,)^)} (5.15) 

up to the second order. Eq. (|5.12|) was used for making the expressions simpler. It is noted 
that the situation here is very similar to that in section 2.2. Applying the same recipe as in 



^^In the bosonic case which we are considering here, we can see no evidence of the spontaneous symmetry 
breaking by the calculation based on the 1 / i'-expansion method and by the numerical simulation as in ref. 
p6t . So we can extend the results here to the N = 00 case. We thank J. Nishimura for informing of this 
fact. 
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section 2.2, we end up with the final expression: 



(P(L)) = C^^^iLY) 
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exp 
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X cos 



iV - 1 / AD 
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(n)- 



(iV + l)LV 
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(iV + l)LV 



1/3- 



D - 1 
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(5.16) 



where up to the first order {n 
The overall constants are 



(1) 



-1) 



Af-l 



X exp 



2 ^ 



7 



(2) 



12 

y 



27V-2JY! 



X exp 



1): 

i(JV2-l)D 

l)Z} + i(iV- 

i{iV2-l)D 
^2 



^, and up to the second order {n 
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At present, we have no result of the computer simulation for the behavior in the region 
of L^g large. However, our result can be expected to represent the behavior correctly at 
least qualitatively. It seems interesting to consider our result together with the asymptotic 
behavior of the eigenvalue density p(A) In fact, p(A) can be obtained from {P{L)) via 

the Fourier transformation: 



P(A) 



dL 



— e-'^^(P(L)). 



-oo 27r 



If we tried to evaluate the integral in the r.h.s. for A large by the saddle point approximation 
using the formula of ( |5.16| ), it would turn out that the width of the Gaussian integral around 
the saddle point becomes infinitely large as A ^ oo. It indicates that we should take into 
account the contribution infinitely far away from the saddle point, and thus that in this case 
the saddle point approximation is inadequate. We need the information of (P(L)) for the 
whole range of L to get the asymptotic behavior of p(A). 

Also, we can consider the improvement of the result in the region of smaller L'^g in the 
previous section by using this method similarly as in the case of the space-time extent in 
the previous subsection. Since it is considered not to change the result essentially, we do not 
present here. 



5.3 Wilson Loop 

We analyze the expectation value of the Wilson loop for the case of the length of the loop 
L large in the improved mean field approximation. Doing the same procedure as before for 
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the operator 



O 



r3i)E^tr( 



we have the following equation for a: 



a^^ D-1 



IM 



(5.17) 



(5.1^ 



{WiL))M 

where the r.h.s. can be easily evaluated by using the expressions ( [4. 281 ) and ( |4.32| ) replaced 
a with a. When L^a^ and L^g are large, the r.h.s. is expanded as 



{WiL))M 2 AT 



N-1 /L^a^^'"^-' ^ 

-e 2 



1 + e 



{my ' V 2 , 

So, the solution up to the next-to-leading order becomes 
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2~2\2Af-3 



(5.19) 
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{N + 1)\{N - 1)\ D 
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+ 0(e-^''^'(LV2)4^-4) , (5.20) 



with cr^ 
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Note that in the L ^ oo limit cr^ coincides to the value of cx^ in the 



N(D-l) ■ 

Gaussian approximation. Up to the first order with respect to S* — S*^/, A is determined as 



1 ^ {{S~S'j^)W{L))m ^ ^..^ 

A (H^(i:))M ' 
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1 + 



A^!(Ar-2)! [ m{N + 3)D N^{N + 2){N + 3) D - 1 
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+0(e-^(LV)2^-^ 



(5.21) 



Plugging these into the unnormalized expectation value {W{L))' = {a'^)2^ e~^' {W{L))m, 
and dividing it by the partition function up to the first order, we end up with the final ex- 
pression: 



{W{L)) 



N+1 



exp 



m-1 
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1 



N^{N + 3)D- 1_ 

4 [ 2 {N + l){N + 6) 1 
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A^!(iV-2)!iV2(A^ + 3) [D N + 2 

+0 [e-'^iL'ay-A + 0{{S - S',,f 



D - 1 



e 2 



(5.22) 
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There are two comments in order. First, in the above expression, {W{L)) has a nonzero hmit 
as L — > oo, which is j^-^exp jj^^j^^j^^i ■ The value approaches to zero as increasing A^. 

Second, since — ^"^^^-^^^^j^ < for > 1, {W{L)) approaches to the hmit from the 
below. In the analysis of the Wilson loop, because a reduces to a in the L — oo limit, the 
result in the Gaussian approximation is considered to roughly reproduce the nonzero limit 
here. In fact, in the L oo limit the result of the Gaussian approximation becomes 



lim 

L— >oo 



(;.'(L))„ - as - S„),y(L))c,„] = ^ + ^,J^^y (5.23) 



which reproduces the expansion of the exponential up to the first order in the limit of (|5.22|) . 
As is seen from the above formula, this seemingly strange behavior becomes hardly visible 
as increasing A^. The phenomenon just like this is reported in ref. DTSf, where for N = 2, 4 



and 8 cases in the supersymmetric IKKT integral numerical simulations are performed and 
results similar as ours are obtained (in Fig. 3a in ref. There, it is interpreted as 

a finite A^ artifact, because of the value of the limit decreasing as increasing A^, which is 
common to our result. Also, its dependence on the shape of the contour is reported in ref. 
T5[] . That is, when changing the contour from the square to a regular polygon, the value 



decreases as increasing the number of edges of the polygon. 

In this method, we can also consider the improvement of the result presented in the 
region of L^g smaller in the previous section. Due to the same reasoning as in the case of 
the Polyakov loop, we do not give here. 



6 Discussions 

Here, we summarize our results and discuss possible future directions for concluding this 
paper. 

We analyzed the reduced Yang-Mills integrals by using the Gaussian approximation and 
its improved version, after confirming the validity of the schemes in the simple example of the 
(/(^-integral. The free energy was evaluated approximately. Since the Gaussian approximation 
can be regarded as the mean field approximation, the approximation is expected to become 
better as the space-time dimensionality D increases. Actually, for the SU{2) case, where 
the exact result is known, we can make sure this expectation. Also for several lower A^'s, 
comparing with the numerical result, we found a tendency of better agreement not only for 
larger D but also for larger A^. Next, we investigated the square of the space-time extent 
and the one-point functions for the Polyakov loop and the Wilson loop both by the Gaussian 
approximation and by the improved mean field approximation. In the SU{2) case, we saw 
that our results for the space-time extent exhibit good agreement with the exact result when 
D is large. In particular, in the improved mean field approximation, the situation becomes 

^^According to M. Staudacher, the authors of ref. have the resuh of the simulation for the Wilson 
loop also in the bosonic case. It shows the similar behavior as in the supersymmetric case ||25|] , which match 
with our result. 
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better than that in the Gaussian approximation. We saw that it hkely holds for general 
by comparing with the 1/D-expansion result. We evaluated the Polyakov loop and the 
Wilson loop by the Gaussian approximation when the length of the loop L is smaller, and 
by the improved mean field approximation when L large. In the former analysis, we saw 



that our results nicely fit the numerical results in ref. [|T^. Furthermore, the remarkably 
simple formulas were obtained in the 't Hooft like large N limit. We also observed quite 
good scaling behaviors in the region > 48, which means that in the scaling region the 
simple formulas represent sufficiently well the behaviors in the case L smaller. In the latter 



analysis, the result of the Wilson loop is in conformity with the numerical result |25|, while 



with respect to the Polyakov loop we do not have any results to be compared with ours as far 
as we know. From the analysis for the simple example, however, our results can be expected 
to reproduce the correct behaviors at least qualitatively. 

For possible future directions, we mainly point out the following two issues. One is an 
extension of this method for supersymmetric systems. Because our motivation is to explore 
approximation schemes trusted nonperturbatively in the IKKT integral, it is one of the most 
important issues. As the first step, it will be good to consider D = 4 supersymmetric case. 



Since in this case there are nice numerical results in ref. |T3], it will be a good test for the 
approximation. From the analysis for various supersymmetric quantum mechanical systems 
in the Gaussian approximation in refs. p2| , it seems to be necessary to consider the 



Gaussian approximation after rewriting the system in terms of the superfields. The other is 
to calculate more general correlators in our framework. In order to use our improved mean 
field formalism for the one-point functions, we considered the restricted class of the operators 
with the isotropic nature. Due to this nature, we could use the simple mean field action. For 
considering general operators, we will need to start with the mean field action with more 
general form. We hope that the analysis along this line gives any insights for studies of the 
nonperturbative aspects of string theory. 
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A Improved Mean Field Analysis in (/)^-integral 

Here we present the improved mean field analysis for 0^-integral as another example which 
holds the assumption mentioned in section 2.2. We consider the 0^-integral defined by the 
classical action: S = 

A.l Saddle Point Approximation 

First, let us evaluate the integral dcj) e~^e^^'^ in the case of L large by the saddle point 
approximation. The solutions for the saddle point equation and the corresponding values of 
—S + iL(j) are 

where k = 1, ■ ■ ■ , 5. Taking into account regions where the integrand does not blow up as 
101 — *• oo, it turns out that we should deform the integration contour so as to pass the three 
points 01, 02 and 03 along each steepest descent direction. After the Gaussian integrals, 
dividing by the partition function 



Z = 

the expectation value becomes 



= liQg'Y^'r (1) , (A.2) 



^ ' ' V 5 6V6r(|)(L(7i/3)2/5 



^2e-|(-¥)(^5^^^)^/^os 



|(sin|)(V/3)6/5_-^ 



(A.3) 



A.2 Improved Mean Field Approximation 



Next, we consider the mean field treatment for this system. For the partition function, we 
start with the following mean field action: 

2 1 /3 

The self-consistency condition (0"^) = (0^)m determines a as = {^j^ . Also, A is given 
by the same form as (|2.18| ). Calculation of the first few terms in the mean field expansion 



leads the expression of the free energy: 

= - In ^ = In/ - ^ In - i - 1^ + 0(((5 - S',,f)c,M), (A.5) 
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where the third and the fourth term in the r.h.s. represent the contribution from {S — S'j^j) m 
and from {{S — S'j^.jY)c,Mi respectively. 

Now let us discuss the improved version of the mean field approximation for the unnor- 
malized expectation value (e*^*^)'. Starting with the mean field action 



2(^2 



1 

+ - 

A 



for the case of O = e*^*^, the self-consistency condition {(f)'^)^ 

. 1 



reads 



X 



(A.6) 



(A.7) 



where we put x = L'^a'^. Since we are interested in the case that x and L^g"^ are large, 
this equation can be solved iteratively. Up to the next-to-leading order, we have the five 
solutions: 



Xk 



Wg' 



+ ^ + 0((LV)-^/^) 



5). (A.8) 



A is given by eq. ( |2.22| ) corresponding to each solution of x. After some calculation with 
respect to the first few terms, we obtain the expression of A as 



A ~ 15 
up to the first order, and as 



+ o{x-') + o{{s-s'^,y) 



1 59 11 ^, 

= — X H h 0(x , 

A 75 25 ^ ' 



OiiS-S[ 



(A.9) 



(A.IO) 



up to the second order. 

The unnormalized expectation value is expressed by the same form as ( |2.27| ) also in this 
case. Now, we compare the exponentiated term — |a; — i in eq. ( p.27|) at each solution {x, A) 
with eq. ( |A.1|) in the saddle point calculation. Then, the following precise correspondence 
similar as in the 0^ case can be found: 

{xi,Xi) ^ (j)2, (a;2, A2) 03, (a;3,A3)^04, (x4,A4)^05, (2:5, As)^^!. (A.ll) 

Thus, we proceed along the same strategy as in the (p'^ case. The solutions (0:3, A3) and 
(x4, A4) lead an unphysical result blowing up as L — > 00, so we discard them. Let us 
combine the contribution from the other solutions with the weight unity. Dividing by the 
partition function, we eventually obtain 

151/6 1 



51/10 (L^l/3)2/5 



^„.5-i/5(sin^)(L(7i/3)6/5_![ 



x[l + Oi{Lg'/Y'/')+0{{S-S',,r^% 



(A.12) 



33 



where up to the first order (n = 1): Ui = t>i = and up to the second order (n = 2): 



_ 3032 „, _ 193 
^2 - 1125' ^2 - ISO- 



Let us compare these with the saddle point resuh ( [A. 31 ). The constant factor in front of 



the whole expression does not exhibit a good result as long as looking at the first two orders: 

||^e~% = 0.36680 ■ ■ ■ up to the 1st order 

ig^e-fri = 0.09029 • • • up to the 2nd order (A. 13) 

>- 



2^ 3 1 _ 0.44819 . . . saddle point method. 



However, the coefficient of the exponential decay or the oscillation, which plays more im- 
portant role with respect to the qualitative behavior, approaches closer to the saddle point 
result as increasing the precision from the first order to the second order: 



29 
30 



. 5-1/5 = 0.70062 ■ ■ • up to the 1st order 



1|| . 5"i/5 = 0.93254 ■ ■ ■ up to the 2nd order (A. 14) 



150 
5 _ 



g — 0.83333 ■ ■ ■ saddle point method. 

B Large Limit in Polyakov loop and Wilson Loop 

Here, we give some detailed explanation about the derivation of ( [4.22|) and (|4.38| ). 



B.l Polyakov Loop 

First, let us start with the Polyakov loop. Eq. ([4.17|) can be rewritten as 



(P(L))o = eW^F(iV + l,2;-^ 



oo 



^--E^(l4)(54)••■(^^)(-)^ 

with a = — . We consider the N ^ oo limit keeping G = y/N q fixed. Because the 

^2{D-l) ^ ^ V y 

series in the r.h.s. converges uniformly with respect to A^, the order of the limit and the 
summation Zl^o t)e changed. So we obtain 

°° 1 1 
hm (P(L))o = E „ ^i„ (-^r = ^^i(2v^), (B.2) 



which is the first term in ( 4.22|) . 



In the second order result ( [4. 211 ), after picking up 0(A^°) -terms dominant in the large A^ 
limit, we have 

]^{{S-S,fP{L))c,o = e^^ 
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F[N + 2,3;- j^]\+0{N-^). 



(B.3) 



Considering the limit of each confluent hypergeometric function as above, we arrive at the 
formula 

jhn ^((^ - SorP{L))c,o = ^p^[^4(2y^) - 3M2V^)]. (B.4) 
This gives the second term in ( [4.22[ ). 



B.2 Wilson Loop 

Next, let us argue about the Wilson loop. We start with considering the limit with respect 
to the function /(L^a^) appearing in the expression of the leading order term ( |4.28| ): 
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We put n = Nx and m = Ny. In the N oo limit, x and y become running over the 
interval 0<x<l,0<?/<l continuously. Then, the first term is expressed as 



/ dx [ dy Jo(2v/xa) Jo(2v^) = \{ I dx Jq{2 
Jo Jo 2 \Jo 



2a 



'M2^)f. (B.6) 



It is seen that the second term can be neglected in the limit, from the following consideration. 
First, we separate the second term into the two parts. One is the contribution from the region 
n — m = 0{N°') (0 < a < 1), and the other is that from n — m = c = 0{N^). With respect 
to the first part, the case of n = 0{N) is dominant, and then ^ < ( const )A^''^", where we 
wrote n — m = hN"- with h being a 0(l)-constant. Also 



m\{{n — m)\] 



< (const) 



e^a \ 1 



2a 



which is negligible in the limit N ^ oo for a fixed arbitrary a. Thus we can neglect the 
first part. Next, let us consider the second part. Noting that n and m run satisfying 
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n — m = c = 0{N^), powercounting leads |(the second part)| < (const)-^ ■ N. So we can 
also neglect this contribution when N oo. 

Thus, the first term alone survives in the limit, and it yields the result 



hm -^JiL'a') = ^iM2y^)y. 
As the consequence, we obtain the limit of the leading order term 



lim {W{L))o = 1 - 



-(Ji(2v/^))2-ll 
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First Order Term For considering the large limit of the first order term, it is convenient 
to rewrite the r.h.s. of eq. ( [4.30|) in terms of connected correlators. Namely, 
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Each connected correlator is given as follows: 
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where 
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Now, let us consider the hmit of the basic quantities fi{L), fsi^L) and iV^/2(L). As the 
result of the argument similar as in the leading order case, for fi{L) we find that only the 
second term in eq. ( [4.35 ) is relevant, and that it gives the limit 



lim ML) = —Ji{2^) 



2^a 



:B.13) 



Also, for fsiL), the term containing the summation of the product of 2nF(—n + 1,2;^°" 
and 2mF{—m, 1; ^^y-) in eq. ( |4.37|) alone becomes relevant. The result is 



2 



dx^/xJl{2^/a) 



L2 



[B.U) 



For N'^f2{L), we go along the same line. With respect to the second, third and last terms 
in eq. (|4.36|) , it can be done easily: 

4a r. 



lim N (the 2nd and 3rd terms in eq. (4.36)) 



N-*oo 



3L2 



{M2y^))^ + {M2V^)f 



lim A^2(the last term in eq. (OBI)) = — r \{Ji{2y/^)f + {J2{2y/^)f 

N^OD 3L2 L 



(B.15) 



In the fourth and fifth terms in eq. ( |4.36| ), the dominant contribution comes from the case 
oin-m = c = 0{N°). 

lim A^^(the fourth term in eq. (4.36)) 



= -ttXI / dx x{Jc{2^/ixx)f + —T-^ i dx^/xJc{2^/xa)Jc+i(2y 
+T^E(c'-c) / dx{U2V^)f 

= ^ [Aa{Jo{2^)f + 4(a - 2){J^{2^)f + A^U2^)J^{2^)] , (B.16) 

where in the last step we used the identity Z^nezl-^nl-^))^ = 1 and some recursion relations 
among Bessel functions. Similarly, 

lim A^^(the fifth term in eq. (4.36)) 



— — ^ / dx xJc^l{2^/axi)Jc+l{2^/xa) = / dx x{Ji{2y/xa)) 

L/ ^_j^Jo L Jo 

-:^,[{M2Mr + {M2V^)f 



:b.i7) 
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Putting eqs. (|B.15|) , ( p.l6|) and ( p.l7|) together, we obtain the simple formula 



hm N'hiL) : 



-Ji(2v^)J2(2yS). 



(B.18) 



Now we can write down the three connected correlators appearing in the r.h.s. of eq. 
(il). The last equation in (|BloD seems to be complicated. However, since as the result of 
powercounting it is only the last term in this equation that gives the leading contribution in 
the large limit, the expression becomes considerably simple. The result is as follows: 

hm N{tTiX,)'W{L))c,o = ^Ji{2V^)J2{2V^) \-{Ji{^V^))' - 1 



lim Ar(tr(XiXiX2X2)W^(L)) 



AT— »oo 



C,0 



3/Vi(2v/^)J2(2 



a 



36 



— {.h{2^))\.h{2^)f + — v^(Ji(2v/^))V2(2v^), 

32 

lim iV(tr(XiX2XiX2)iy(L))c,o = 77(^2(2v^))^ 

Plugging these into ( |B.9|) , we end up with the final expression: 

1 2 



(B.19) 



-{iS-So)WiL))c,o 



9(Ji(2v/^))2(J3(2v^))2 + 8(J2(2v^))^ 



Eqs. ( |B.8| ) and ( |B.20| ) give the r.h.s. of eq. ( [4.38| ). Thus we complete the derivation. 



B.20) 
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